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ILLUSTRATIONS 
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I INTRODUCTION 


The  theoretical  work  presented  in  this  report  was  undertaken  with 
the  long-range  objective  of  developing  a better  basic  understanding  of 
the  initiation  and  propagation  of  detonation  in  condensed  explosives. 

The  achievement  of  this  objective  is  necessary  in  a program  for  the 
improved  characterization,  control,  and  effective  use  of  explosives  in 
military  applications.  Theoretical  studies  are  an  essential  part  of  such 
a program  because  experimental  studies  of  detonating  explosives  are  not 
completely  definitive  and  depend  on  models  of  shock-induced  reactive  flow 
for  their  analysis  and  interpretation. 

1-3 

Studies  of  Zeldovich-von  Neumann-Doering  (Z-N-D)  waves  were 
performed  to  elucidate  and  determine  the  parameters  that  govern  the  ini- 
tiation of  detonation.  The  problem  of  a single  shock  trajectory  for 
buildup  to  detonation  was  considered,  and  the  differential  equation 
governing  a shock  discontinuity  was  examined  to  determine  different  con- 
ditions for  this  type  of  flow.  One  of  these  conditions  was  used  to  model 
the  type  of  flow  observed  in  PBX  9404  during  the  early  stages  of  initiation 
produced  by  a flying  plate. 


Other  properties  of  initiation  produced  by  a constant  velocity  piston 
were  considered.  The  equations  relating  the  initial  flow  to  the  initial 
energy  release  rate  in  such  a wave  were  derived,  as  well  as  conditions 
for  the  shock  to  accelerate  with  either  a positive  or  a negative  pressure 
gradient.  These  conditions  demonstrate  the  dependence  of  the  initiation 
process  on  the  energy  release  rate,  the  sound  speed,  and  the  relationship 
between  these  quantities.  A critical  energy  was  defined  for  waves  that 
build  up  to  detonation  with  a positive  pressure  gradient. 

Work  on  reactive  shocks  was  performed  with  the  objective  of  determin- 
ing the  dependence  of  initiation  on  the  energy  release  rate  and  the  rear- 
boundary conditions.  Exact  integral  relationships  for  unsteady  flow  were 
derived  as  generalized  Rankine-Hugoniot  equations,  but  the  attempts  to 
construct  explicit  solutions  for  the  buildup  to  steady  detonation  in  a 
polytropic  material  with  a prescribed  energy  release  rate  were  unsuccessful. 
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II  BASIC  ASSUMPTIONS 


Although  the  basic  assumptions  of  our  theoretical  studies  are  well 
defined  in  reference  4,  they  will  be  repeated  here  for  completeness. 
Reactive  shocks  are  treated  as  (Z-N-D)  waves  with  no  reaction  across  the 
shock  discontinuity,  and  the  flow  induced  by  the  shock  is  assumed  to  be 
adiabatic,  inviscid,  and  in  local  equilibrium.  Consequently,  irrevers- 
ible processes  are  restricted  to  shock  compression  and  chemical  reactions 
behind  the  shock. 

Let  e,  p,  and  P = 1/v  denote,  respectively,  specific  internal 
energy,  pressure,  and  density--the  inverse  of  the  specific  volume.  Then 
the  Rankine-Hugoniot  jump  conditions  governing  the  shock  discontinuity 
propagating  into  stationary  material  at  the  front  of  the  wave  can  be 
written  as: 


PH(U  - V 

= PoU 

(i) 

PH  - P° 

■ P°UUH 

(2) 

eH-e° 

= £(P„  + PoHvo  - v ) 

H tl 

(3) 

where  U denotes  shock  velocity,  u denotes  particle  velocity,  and  the 
subscripts  H and  o denote  quantities  immediately  behind  and  in  front 
of  the  wave.  Since  our  consideration  will  be  restricted  to  rectilinear 
flow,  it  is  convenient  to  write  the  conservation  equations  governing  the 
flow  behind  the  shock  as: 


1 dv 

du 

(4) 

v dt 

dr 

h 

t 

du 

dp 

(5) 

Po^ 

""  Sh 

h 

t 

de 

dv 

(6) 
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where  h and  r denote  Lagrange  and  Euler ian  distance  coordinates  and 
t denotes  the  time.  The  flow  and  the  chemistry  are  coupled  by  the 
dependence  of  e on  the  reaction  coordinates.  For  convenience  in  this 
report,  only  a single  exothermic  reaction  will  be  considered. 


Let  s,  T,  A,  and  X denote  the  specific  entropy,  temperature,  the 
chemical  affinity,6  and  the  extent  of  reaction.  Then  Gibbs'  equation  can 
be  written  as 


(7) 


and  the  combination  of  Eqs.  (6)  and  (7)  gives  the  corresponding  equation 
for  the  entropy  production  in  shocked  material  as 


9s  dX 

at  = A at 

h 


Combination  of  Eq.  (8)  with  the  time  derivative  of  the  s = s(p,  P,  X) 
equation  of  state  leads  to  the  equation 


dp  a dp  r dX 

dt  C dt  + v q dt 

h h h 


(9) 


where  the  frozen  sound  speed  c = (dp/dp) s ^ , the  Griineisen  parameter 
T/v  = (dp/de)Vf^,  and  the  specific  heat  of  reaction  q = -(de/dX)p,v. 
Equation  (9)  is  the  most  convenient  equation  for  studying  the  interaction 
between  the  chemistry  and  the  flow  in  the  present  treatment  of  Z-N-D  waves. 
It  can  be  derived  more  directly  by  combining  Eq.  (6)  with  the  time 
derivative  of  the  e = e(p,  v,  X)  equation  of  state. 
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Ill  EQUATIONS  OF  STATE 


Equations  of  State  of  Explosives 


Equations  of  state  of  condensed  explosives  can  be  constructed  from 

the  results  of  shock  wave  experiments  that  have  been  performed  at  the 

T 8 

Ballistic  Research  Laboratories  and  the  Naval  Ordnance  Laboratory  over 
the  last  few  years.  These  investigations  lead  to  the  conclusion  that  the 
unreacted  Hugoniot  curve  of  a condensed  explosive  can  be  adequately 
described  by  a linear  relationship  of  the  form 

U = a + bu  (10) 

H 

in  the  (U-u  ) plane.  The  values  of  the  constants  a and  b for 
different  explosives  are  given  in  References  7 and  8.  Combination  of 
Eq.  (10)  with  Eqs.  (1)  and  (2)  gives  the  corresponding  Hugoniot  in  the 
(p-v)  plane  as 

a2(v0  - v) 

p = (11) 

[v0  - b (vo  - v)] 


and  the  variation  of  specific  energy  along  the  Hugoniot  is  determined 
by  Eqs.  (11)  and  (3).  Differentiation  of  Eq.  (11)  gives  the  slope  of 
the  Hugoniot  in  the  (p-v)  plane  as 


dp  _ a [v0  + b(v0  - v)3  . 

dv  a [v0  - b (v0  - v)]3 

H 

and  the  corresponding  equation  for  the  slope  of  the  Hugoniot  in  the 
(p-v)  plane  is  obtained  from  the  jump  conditions  as 

= P0(a  + 2 b u)  (13) 

du 

H 
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2 

An  expression  for  the  sound  speed  along  the  Hugoniot  (c/v)  in 

terms  of  (dp/dv)  can  be  obtained  by  combining  the  Hugoniot  equation 
n 

with  the  first  and  second  laws  of  thermodynamics.  Thus  combination  of 
the  equations 


(14) 


- £p  + £(vo  - v) 


dp  _ /c\2  Ft  ds 

dv  \v / v dv 

H H H 


(15) 


(16) 


to  eliminate 


Tds/dv 

H 


leads  to  the  equation 


dp  (T  / 

dv  L2  W 
H 


(17) 


An  assumption  about  T allows  the  (e-p-v)  equation  of  state  of 
an  explosive  to  be  constructed  over  the  volume  spanned  by  the  Hugoniot 
curve.  If  r is  assumed  to  be  constant,  or  a function  of  v,  for 
example,  the  (e-p-v)  equation  of  state  has  the  form 


e = p-  + g(v) 


(18) 


which  can  be  rewritten  as 

* - CH  = f (P  ’ P«)  <19> 

since  e„  ,j:d  p„  are  known  as  functions  of  v along  the  Hugoniot 
ft  H 

curve.  An  additional  assumption  about  the  specific  heat  at  constant 

volume  C that  is  consistent  with  [9(F/v)/9p]  =0  is  then  required 

v v 

to  construct  the  corresponding  (p-T-v)  equation  of  state. 


Equations  of  State  of  the  Explosive  Mixture 


The  equation  of  state  of  the  explosive  mixture  will  be  considered 
in  more  detail  because  it  may  be  necessary  to  remove  the  assumption  of 
local  thermal  equilibrium  at  a later  date. 


Let  the  subscripts  2 and  1 denote  the  explosive  and  its  products, 
and  let  the  superscript  0 denote  the  standard  state.  Then  the  equations 
for  the  specific  energy,  specific  entropy,  and  specific  volume  of  the 
explosive  mixture  can  be  written  as 


e|  - Aq  + Aei  + (1  - A)ea 

(20) 

Asx  + (1  - X)S2 

(21) 

Avi  + (1  - A)va 

(22) 

where  ei  denotes  the  heat  of  formation  of  the  explosive,  and  Q = ei  - e? 
denotes  the  standard  heat  of  reaction.  If  ex  and  e2  are  considered  to 
be  functions  of  mechanical  variables,  then  these  equations  of  state  are 
written  formally  as  ei  = ei(p,vi)  and  es  = es(p,v2)  with  pi  = ps  = p 
because  the  mixture  is  assumed  to  be  in  mechanical  equilibrium.  But  if 
pressure  and  temperature  are  chosen  as  independent  variables,  they  will 
be  written  as  ei  = ei(p,Ti)  and  es  - ea(p,Ta).  The  assumption  of 
local  thermal  equilibrium  is  then  expressed  by  the  condition  Ti  = Ta. 

Note  that  Eq.  (10)  simplifies  to  the  equation  given  in  reference  4, 

e = ei  - AQ  + (23) 

when  the  explosive  and  its  products  are  assumed  to  have  the  polytropic 
equation  of  state  with  the  same  value  of  the  polytropic  index  K. 


Combination  of  Eq,  (6),  expressing  the  First  Law  of  Thermodynamics, 
with  the  time  derivative  of  Eq.  (20)  leads  to  the  equation 


Bex 

3t 


+ (1  - X) 


Be  2 
3t 


[Q  + (e2  - ex)]  ~ - 


Bv 

p bT 


(24) 


where  the  Lagrange  subscripts  have  been  omitted  for  notational  simplicity. 
The  ex(p,vx)  and  ea(p,va)  equations  of  state  can  then  be  used  to 
transform  the  left-hand  side  of  Eq.  (24)  into  the  expression 


♦ ] * W&)  It 

v a P 


(25) 
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Setting  X = 0,  and  dv/dt  = dv2/dt  + (v2  - vi)dX/dt  in  Eq.  (24),  with 
the  left-hand  side  written  as  shown  in  (25) , gives  the  following  equation 
for  conditions  at  the  shock  front 


/deaN 

dp 

r r 

Up  / 

at  + 

A aw  + p3ti 

r = L( 

(2 


It  follows  from  Eq.  (26)  that  the  values  of  ei  and  v*  are  in  general 
required  to  calculate  conditions  at  the  shock  front  even  though  the 
reaction  is  just  starting.  For  the  mixture  governed  by  Eq.  (23),  how- 
ever, only  Q appears  on  the  right-hand  side  of  Eq.  (26)  because 
ei  - pv i / ( K - 1),  e2  = pv2/( t£  - 1),  and  e2  - ei  = (dea/dv2)(V2  - Vi). 


When  Ti  = T2,  vi  and  V2  can  be  considered  as  functions  of  v, 
p,  and  X because  the  T = T(p,v,X)  equation  of  state  of  the  mixture 
can  be  formally  substituted  into  the  vi(T,p)  and  V2(T,p)  equations 
of  state  of  its  constituents.  In  this  case, 


avi 

/9v i 

\ 

/^viN 

ap 

/ avi 

\ ax 

at  = 

\dv 

/ .at 

P/ 

Up  ) 

V, 

\ax 

/ at 

P.V 

aV2 

/5va 

/dva\ 

Ip 

/dva 

\ ax 

at  " 

\bv 

).bt 

Up  / 

. at 

\ax 

/ at 

rX 

v, 

X 

P»  v 

and  Eq.  (24)  can  be  written  as 


(27) 


(28) 


with 


- *(£)  * (It)  (ir)  J + (l  - »[(£ ) * (IfXIr)  J 

Vl  p V,X  V2  V,  X 


Differentiation  of  Eq.  (22)  gives  the  relationship  for  the  volume  deriva- 
tives in  Eq.  (32)  as 

) + (1  - *>(f^)  “ “<vi  " Va)  <33) 

p,v  P,v 


and  Eq.  (32)  reduces  to  the  equation  (3e/3X)  = - Q for  the  mixture 

satisfying  Eq.  (23). 
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IV  A SINGLE  SHOCK  TRAJECTORY 
I OR  BUILDUP  TO  DETONATION 


The  use  of  a single  shock  trajectory  to  describe  the  initiation 

g 

of  detonation  in  a condensed  explosive  was  postulated  by  Lindstrom. 
Lindstrom' s postulate  is  considered  here  because  the  existence  of  such 
a single  shock  trajectory  leads  to  interesting  conclusions  about  the 
initiation  process.  Our  consideration  is  based  on  the  following  differ- 
ential equation4  governing  the  pressure  along  the  shock  path 


where  (du/dp)^  denotes  the  slope  of  the  nonreactive  Hugoniot  curve  in 
the  (u-p)  plane.  It  is  convenient  to  rewrite  Eq.  (34)  for  a strong 
shock  in  material  governed  by  Eq.  (23).  In  this  case,  the  mechanical 
jump  conditions  can  be  written  as 


t 


U 


K + 1 
2 


u 


H 


(35) 


(38)  I 

(37) 

I 


(34)  simplifies  to  the  equation 


where  subscript  D denotes  the  state  at  the  front  of  the  steady-state 
Chapman-Jouguet  (CJ)  wave,  T = (K  - 1),  and  q and  p are  related 
by  the  expression  (K-l)Pq=(K+l)p/4(K-l),  J 


I 


Equation  (38)  cannot  in  general  be  integrated  because  the  reactive 
flow  is  subsonic  with  respect  to  the  shock,  and  the  particle  velocity 
gradient  at  the  front  of  the  wave  depends  on  the  flow  between  the  shock 
and  the  rear  boundary.  There  will  be  a unique  shock  trajectory  in  either 
of  two  cases.  In  the  first,  the  particle  velocity  gradient  at  the  shock 
front  depends  on  a shock  parameter.  In  the  second,  the  piston  motion, 
initiating  the  wave  with  initial  particle  velocity  u < u^  , is  itself 
generated  during  the  buildup  phase  of  the  wave  Initiated  with  a lower 
initial  particle  velocity  u^ / < u^.  The  first  case  will  be  discussed 
before  the  second. 

When  the  particle  velocity  gradient  at  the  front  of  the  wave  depends 
on  a shock  parameter,  Eq.  (38)  becomes  an  ordinary  differential  equation. 
It  can  then  be  integrated  formally  to  give  the  following  type  of  equation 
for  the  dependence  of  shock  pressure  on  time, 


p = p P(t/Gf)  (39) 

H C 

where  p denotes  the  minimum  initial  pressure  for  the  onset  of  detona- 
tion, pfo)  = 1,  and  Of  is  a characteristic  reaction  time.  The  intro- 
duction of  a critical  pressure  for  the  onset  of  detonation  is  reasonable 
because  many  explosives  are  found  to  be  nonreactlve  below  about  40  kbar. 
For  waves  that  build  up  to  detonation,  the  function  P(t/0f)  also  satis- 
fies the  condition  P(t/tt)  p^/p  as  t /ft  "*  °°.  The  relationship  for 
shock  pressure  in  a wave  with  an  initial  shock  pressure  p^  > p^  can  be 
written  as 

PH  = PiPt(t  " T)/tt]  (40) 

with 

P.  = P P(T/0)  (41) 

1 c 

A simple  solution  to  Eq.  (38)  is  constructed  to  exemplify  this  type  of 
behavior.  The  reaction  rate  along  the  shock  path  is  assumed  to  satisfy 
the  conditions 
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It  is  convenient  to  introduce  the  dimensionless  variables  P = p^/p^ 

and  P « p /p  , and  rewrite  the  differential  equation  as 
c c D 


dP  (K  + 1) 

r.p  - p . 
( c> 

1 r (**)  1 

dt  3 

u - p > 

L c 

2(K  - 1)0^  Vdr^J 

(43) 


l 


I 

i 

; 


It  is  now  necessary  to  prescribe  (du/Sr)^  as  a function  of  the 
shocked  state  to  obtain  a single  shock  trajectory.  The  gradient  should 
also  be  prescribed  so  that  it  attains  the  correct  value  in  the  CJ  wave. 
For  simplicity,  the  particle  velocity  gradient  was  assumed  to  be 


du 

dr 

H 


, (P  - P ) 
1 c 


2(*  - 1)0^  (1  - Pc) 


(44) 


so  that  (du/dr)  ~*  (du/dr)^  as  P “*  1,  and  the  shock  propagates  at 
constant  velocity  when  either  P = P^  or  P^  = 1.  The  combination  of 
Eqs.  (43)  and  (44)  gives  the  differential  equation 


(1  - P ) 
c 


dP 


(K  + 1) 


(P  - P )(1  - P)  dt  6(K  - 1)0^ 
which  can  be  integrated  to  give 


(P  - P )(1  - P ) = (P  - P )(1  - P)  exp  [(K  + l)t/6(K 

c i i c 


(45) 


- 


(46) 


Equation  (46)  with  K = 3 was  used  to  plot  the  graph  of  time  to  detona- 
tion versus  initial  pressure  P shown  in  Figure  1.  The  dimensionless 
time  to  detonation  t*/0^  was  defined  as  the  time  taken  for  P to 
attain  a value  of  0.95. 


Consider  now  the  case  when  a piston  initiating  the  flow  at  p > p 
and  a particle  with  the  same  initial  pressure  in  the  wave  initiated  at 
p follow  the  same  path.  Such  a situation  will  occur,  for  example, 
wfien  the  acceleration  of  a particle  depends  only  on  its  initial  pressure 
or  particle  velocity.  It  will  also  occur  when  the  acceleration  is  zero 
and  the  particle  paths  are  straight  lines.  An  exact  solution  for  this 
type  of  flow  is  presented  in  this  report  because  it  was  observed  by 
Kennedy10  in  the  early  stages  of  initiation  of  detonation  in  PBX  9404. 
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FIGURE  1 PLOT  OF  TIME  TO  DETONATION  VERSUS  INITIAL  PRESSURE 


In  these  experiments,  flying  plates  thrown  at  different  velocities  behave 
as  constant  velocity  pistons  and  produce  the  same  shock  trajectory  when 
the  energy  release  rate  is  small.  To  satisfy  the  momentum  equation, 

Eq.  (5),  the  pressure  must  be  a function  only  of  the  time  t , and  the 
particle  velocity  must  be  a function  only  of  the  Lagrange  distance  h. 

It  is  convenient  to  introduce  the  Lagjrange  time  coordinate  T by  the 
equation  for  the  shock  path  h(T)  = J Udt , so  that  T is  the  time  a 
particle  enters  the  wave  and  t = T along  the  shock  path.  The  shock 
pressure  p and  the  shock  particle  velocity  uR  can  then  be  used  to 
define  the  rime  dependence  of  the  pressure  and  particle  velocity  fields 
ns  P - p^(t)  and  u = u^ ( t ) , and  the  Lagrange  continuity  equation, 
written  as 

p(l?)  - VB(T>  - YT)I  (47> 

can  be  used  to  obtain  the  corresponding  density  field. 

The  poly tropic  equation  Eq.  (23)  will  again  be  used  for  convenience. 
In  this  case  the  shock  is  governed  by  the  equations 


(48) 

(49) 

(50) 


dp  „ du 

StH  = ' P-arH 


<*  - 1,Vh 


( 


»*  7 ‘ '' 


and  <JH  = Q dX/dtg.  When  the  particle  paths  are  straight  lines  and 
(dp/dh)^  = (du/dt)T  * 0 , combination  of  Eqs.  (48)-(51)  givea  the  rela- 
tionship between  the  particle  velocity  gradient  and  the  energy  release 
rate  at  the  shock  front  as 


and  the  corresponding  shock  trajectory  Is  determined  by  the  equation 


du 


u. 


H 


H dt 


■ <.K  ~ I)  Q 
(2*  - 1) 


Integrating  Eq.  (53)  gives  the  particle  velocity  field  as 
u(T)  = ujl  + Ki  I(T)]* 


(53) 


(54) 


and  the  pressure  field  follows  from  the  jump  conditions  as 

P(t)  = Pitl  + Ki  I(t)]  (55) 

T=t 

where  Ki  = 2(K  - 1)/(2K  - 1),  u®  I(T  = t)  - J Q^s^ds',  and  the 
subscript  i denotes  the  initial  condition  at  0 t » T m 0. 

Differentiating  £q.  (54)  gives  the  Lagrange  particle  velocity  gradient 
as 


du 

dT 


Ki 

2 


fyT) 

U(T) 


(56) 


integrating  Eq.  (56)  gives  the  deformation  gradient  as 


dr 

dT 


(K  - 1)U®(T)  + Kx6  (T)(t  - T) 
n H 

U(T) 


(57) 


and  the  Eulerian  particle  velocity  gradient  follows  from  Eqs.  (56)  and 
(57)  as 


du  ^ K^H<T)  

at  " (K  - l)u*(T)  + K1QH(T)(t  - T) 


(58) 


Equation  (57)  and  the  Lagrange  continuity  equation  give  the  equation  for 
the  specific  volume  as 


v 


KiQH(T)(t  - T) 

1 + <K  - l)u*(T) 
n 


(59) 


ft':STY  &&  I -jfi-  M i ♦.  . A «■ 

It  follows  from  Eq.  (53),  (54),  (55),  and  (59)  that  the  energy  release 
rate  at  the  shock  front  determines  the  flow.  The  case  when  the 


particle  velocity  gradient  is  a function  of  time  is  of  particular  inter- 
est and  will  be  considered  in  more  detail.  Equation  (58)  gives  the 
condition  for  this  type  of  flow  as 

(tC  - l)u*(T)  - KiOgC^HT  + a)  - 0 (60) 

where  a is  a characteristic  tine  defined  by  the  initial  reaction  rate, 
and  the  corresponding  equations  for  the  particle  velocity  gradient  and 
the  shock  particle  velocity  follow  from  Eqs.  (58),  (60),  and  (53)  as 


du  _ 1 

dr  ~ t + ot 

and 

% _ (K  - 1)  UH 
dT  “ 2 (T  + 00 


(61) 


(62) 


Integration  of  Eq.  (62)  gives  the  particle  velocity  field  as 


K-l 

u_  _ + on~ 

u±  L a J 

the  pressure  field  as 


(63) 


and  the  density  field  follows  from  Eq.  (59)  as 


(64) 


P_ 

P. 


H 


- GtS 


The  energy  equation 


dp  du 

3t  " p S + ' 


- 1)P$ 


obtained  by  setting  ca  » Kp/p  and  Ri 

the  volumetric  energy  release  rate  as 

K-B 


PQ 


(2IC  - 1)  h (eJT* 

(1C  - 1)  tt  \p  J 


(65) 


(66) 

(K  - 1)Q  in  Eq.  (9),  gives 

; A wS  i:"s  ■ v 

J$P  ij  -v  if. 

(67) 
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Since  the  shock  pressure  and  particle  velocity  increase  with  time, 
Eqs.  (63) , (65)  and  (67)  cannot  be  used  to  describe  the  whole  of  the  ini- 
tiation process  resulting  in  a steady  detonation  wave.  Other  equations 
must  therefore  be  patched  to  them  to  describe  the  approach  to  the  steady 
state.  It  is  reasonable  to  assume  that  reaction  in  the  early  stages  of 
initiation  occurs  in  hot  spots  arising  from  inhomogeneities.  Then  the 
energy  release  rate  in  hot  spots  produced  at  the  shock  front  determines 
the  flow,  and  the  associated  energy  release  rate  behind  the  accelerating 
shock  will  enhance  reaction  in  the  bulk  of  the  shocked  material.  If, 
moreover,  there  is  a critical  temperature  for  the  onset  of  thermal 
explosion,  the  temperature  distribution  in  the  wave  will  determine  where 
shocked  material  explode*  and  produces  a second  shock,  which  initiates 
bulk  reaction  and  dominates  the  flow  in  the  later  stages  of  the  initiation 
process. 


Consider,  for  example,  the  simple  case  when  (dT/dX)  = (dp/dX)_,  = 0, 

(de/dX)p  y = (de/dX)T  w , and  the  time  rate  of  change  of  temperature  along 
a particle  path  and  tfte  Lagrange  temperature  gradient  are  governed  by  the 
equations. 


dT  _ S/  ^ dp  /dT\  dv 

dt  v dt  + \dv/  N dt 

p,X 


(68) 


and 


dT 

dT 


(69) 


Combining  Eqs.  (68)  and  (69),  subject  to  the  condition  (dv/dt)H  = 
-(dv/d'Oji  along  the  shock  path,  gives  the  equation  governing  the  shock 
temperature  as 


% _ V*  - » % 

dt  v dt 


(70) 


since  dp/dt  * dpjj/dt,  Since  dv/dt  >0  it  follows  from  Eqs.  (68)  and 
(70)  that  the  temperature  increases  faster  along  a particle  path  than 
along  the  shock  path.  The  mechanism  of  initiation  will  resemble  that 
observed  in  nitromethane11  when  the  initial  pressure  is  about  80  kbar, 
because  the  critical  temperature  will  be  attained  and  thermal  explosion 
will  occur  behind  the  first  shock.  The  position  where  the  second  shock 
is  formed,  however,  depends  on  the  temperature  distribution  and  can  be 
determined  when  the  explicit  form  of  the  p ■ p(v,T)  equation  of  state 
is  known. 
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V INITIATION  BY  A CONSTANT  VELOCITY  PISTON 


Initiation  by  a constant  velocity  piston  was  considered  because 
flying  plates  are  used  in  experimental  initiation  studies  and  because 
constant  particle  velocity  is  chosen  as  the  rear  boundary  in  many  numeri- 
cal studies.  The  problem  is  to  understand  and  quantitatively  account  for 
the  observed  shock  trajectories  in  condensed  explosives.  Of  particular 
interest  is  the  case  when  the  initial  shock  pressure  is  low  and  the  shock 
propagates  at  essentially  its  initial  velocity  before  accelerating  very 
rapidly  to  attain  the  detonation  velocity.  Since  the  wave  can  develop 
initially  with  either  a positive  or  a negative  pressure  gradient  but 
must  have  a positive  gradient  in  the  final  buildup  to  detonation,  our 
approach  was  to  determine  a condition  for  the  initial  shape  of  the  wave 
and  then  look  for  a criticality  condition  for  the  onset  of  detonation. 

It  is  convenient  first  to  derive  the  initial  conditions  in  a wave  ini- 
tiated by  a piston  with  an  arbitrary  motion.  The  combination  of  Eqs. 

(48) -(50)  with  Eq.  (9)  evaluated  at  the  shock  front  gives  the  following 
equation  for  the  derivatives  of  the  particle  velocity  along  the  shock 
path: 


dp  1 du 

dV  % 


r ax 


(71) 


It  follows  from  the  momentum  equation,  Eq.  (5),  that  Eq.  (71)  is  also  an 


equation  relating  the  pressure  and  particle  velocity  gradients  at  the 
shock  front.  The  initial  particle  velocity  gradient  in  a wave  initiated 


by  a constant  velocity  piston  is  obtained  by  setting  du/dt  = 0 in 
Eq.  (71)  as 


v[U(dp/du)  + v0(c/v)  ] 


where  q^  ■ q 3X/dt^.  When  Ou/dh)  is  known,  the  initial  time  rate 
of  change  of  pressure  in  the  wave  ana  the  initial  time  rate  of  change 
of  shock  particle  velocity  can  be  evaluated  with  the  equations 


and 


du 

H 

dt 


(74) 


H 


obtained  by  combining  Eqs.  (48)- (50)  with  the  boundary  condition 
du/dt  = dp/dh  ■ 0,  The  corresponding  value  of  the  volume  gradient  can 
be  evaluated  with  the  equation, 


dv  fdv  vq]  du 

dh  “ Ldu  U J dh 

H H H 


(75) 


obtained  by  combining  the  continuity  equation  with  the  identity 


dv  dV>H  _ dv  dv 

du  dt  dt  + dh 

H H H 


(76) 


It  follows  from  Eqs.  (72)-(7S)  that  the  initial  values  of  the  first 
derivatives  of  the  flow  in  a reactive  wave  initiated  by  a constant 
velocity  piston  are  determined  by  the  initial  shock  strength,  the  equa- 
tion of  state  of  the  explosive,  and  the  initial  value  of  the  energy 
release  rate.  Alternatively , a knowledge  of  the  initial  value  of  one 
of  the  derivatives  is  sufficient  to  determine  the  initial  value  of  the 
energy  release  rate  when  a,  b,  and  T are  known.  Measurement  of  one 
of  these  quantities  in  a series  of  experiments  using  flying  plates 
thrown  at  different  velocities  is  therefore  sufficient  to  determine  the 
dependence  of  the  initial  energy  release  rate  on  pressure. 

We  will  now  derive  the  equations  that  determine  the  shape  of  the 
wave  during  the  early  stages  of  the  initiation  process.  For  notational 
convenience,  the  subscript  H will  now  be  omitted  from  the  derivatives, 
but  it  should  be  remembered  that  the  equations  specify  the  initial  con- 
ditions in  the  wave.  Since  the  piston  moves  with  constant  velocity,  it 
follows  from  the  momentum  equation  that  the  pressure  gradient  at  the 
piston  is  zero.  Consequently,  the  initial  rates  of  change  of  pressure 
along  the  piston  and  along  the  shock  path  are  equal  as  shown  by  Eq.  (73), 
and  it  is  necessary  to  find  the  second  derivative  (d  p/dh  ) to  deter- 
mine whether  the  wave  develops  with  a positive  or  a negative  pressure 
gradient.  During  the  early  stages  of  the  flow,  the  wave  will  develop 
with  a positive  gradient  and  the  pressure  will  be  higher  at  the  shock 
than  at  the  piston  when  the  initial  value  of  d8p/dh*  > 0;  when  dap/dh*  < 0, 
the  wave  will  develop  with  a negative  gradient  and  the  pressure  will  be 
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higher  at  the  piston  than  at  the  shock;  when  3 p/3h  = 0,  there  will 
be  no  gradient. 


The  equation  for  (3ap/3ha)  is  derived  from  Eqs.  (4),  (5),  (9), 
and  Eqs.  (48-50).  Combination  of  the  equations  obtained  by  differen- 
tiating Eqs.  (4),  (5),  and  (9)  to  eliminate  the  second  derivatives  of 
shock  pressure  and  particle  velocity  leads  to  the  equation 


U 


( 


U + 2v0 


dp\  _ 3 p dp 
du / 3t^  + du 


3au  3u  dU\  dap  / dUH\2 
3hs  + 3h  dt/  + dua  \dt  / 


(77) 


when  account  is  taken  of  the  equations 


3au  _ 3ap 

3hdt  V°  3ha 


(78) 


*2  >2 

o u o p 

3ta  3t3h 


(79) 


imposed  by  the  rear  boundary  condition.  Partial  dif ferentation  of 
Eq.  (9)  written  in  terms  of  3u/3h,  with  respect  to  t and  h,  yields 
expressions  for  3ap/3ta  and  3au/3ha.  Substitution  of  these  expressions 
into  Eq.  (77)  then  gives  the  required  expression  for  3 p/3h  . 

2 2 

The  equation  obtained  for  3 p/3t  is 


aaP 

3^ 


f cvo\  3 p 

3u  3 /< 

cN3  r/v 

of  Su . r 

( V / Sh*  ” 

V°  3h  3t  v 

;)  - r(.t 

/ V°  3hq  + v q 

(80) 


2 2 

and  the  equation  obtained  for  3 u/3h  is 


3au  _ fvf  3u  3 /c\a  /v\a_jY  vo\  1 3q  1 dv  q *| 
3hy  \c/  3h  3h\v/  + \c/  L\  v / v^  3h  v7  oh  v0J 


(81) 


3 2 

The  combination  of  Eqs.  (70),  (80),  and  (81)  to  eliminate  3 p/dt  and 
3 3 2 3 

3 u/3h  give  the  equation  for  3 o/3h  as 
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- lu:  1 


}h~:  - 


+ 


B = s [T° 


+ 


where 


a dp  (v U f 5 fcf  dp  dU  d*p  dUHl 

V°  du  Vvqc/  Bh  \v/  + du  dt  + du^  dt  J 

Iq  r /Vof  Bu  dp  /u'f  Bv"j 
v0  L\v  / 3h  du  \c/  BhJ 


(82) 


B = [u2  + 2Uv0dp/du  + (vqc/v)8] 


Equation  (82)  specifies  the  dependence  of  the  initial  shape  of  the 
pressure  profile  produced  by  a constant  velocity  piston  on  the  energy 
release  rate  and  the  e = e(p,  v,  X)  equation  of  state,  and  can  be 
used  to  determine  whether  the  wave  develops  with  a negative  or  a positive 
pressure  gradient  when  these  quantities  are  known. 


Consider  the  case  when  the  explosive  is  described  by  Eq.  (18),  the 
products  obey  a polytropic  equation  of  state,  and  Eq.  (20)  gives  the 
e = e(p,  v,  X)  equation  of  state  as 


e = e%  - Aq  + tc”--  + (1  - X) [Kapva  + g(va)]  (83) 


with  Ka  = (I*”1  - (K  - 1)  *).  Note  that  Eq.  (83)  reduces  to  Eq.  (23) 
when  the  explosive  is  assuaed  to  have  a poly tropic  equation  of  state 
since  then  Ka  = g(va)  = 0.  Equation  (82)  can  be  used  to  derive  expres- 
sions  for  the  derivatives  3(c/v)  /dt  and  3(c/v)  /Bh  that  are  needed 
to  determine  the  sign  of  3 p/3h  . The  simplest  case  with  the  explosive 
and  its  products  in  thermal  equilibrium  (Ti  = Ta  = T)  was  considered 
so  that  the  identity 


p + (3e/3v)  . 
PA 

<a*/ap)v,X 


(84) 


3 3 

could  be  used.  It  is  clear  that  the  second  derivatives  (B  e/Bp  ).  , 

(B^/BpBv).^,  (BV/Bv2)^,  (B^/BXBp)^,  and  (B^/BXBv)^ 

must  be  calculated  in  order  to  calculate  the  derivatives  of 
(c/v)8  from  Eq.  (84).  Differentiation  of.  Eq.  (82)  leads  to  the 
following  equations  for  these  derivatives 
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s 0 


(85) 


(86) 

(87) 

(88) 

(89) 


Equations  (88)  and  (89)  define  the  influence  of  the  reaction  on  the 
sound  speed  when  the  explosive  and  its  products  satisfy  Eq.  (83).  Fur- 
thermore, it  follows  that  Eq.  (83)  is  in  general  insufficient  to  calcu- 
late a2p/3h2  when  the  energy  release  rate  is  known.  Values  of 
d3p/dh3  based  on  additional  equation  of  state  assumptions  were  calcu- 
lated and  will  be  discussed  later.  But  most  attention  was  given  to  the 
simple  case  when  the  explosive  is  polytropic  and  the  sound  speed  is  not 
influenced  by  the  -reaction  because  all  the  second  derivatives  of  the 
energy  in  Eqs.  (85)- (89)  vanish  except  3se/5pdv  *=  1/(K  - 1). 

Evaluation  of  the  terms  in  Eqs.  (81)  and  (82)  using  Eq.  (23)  gives 
the  following  equations 


a d2u 
X1  Sh5 


= 2 r (.K..3_i)  j&  1 

It  L(K  + 1)  3h  ill 


2(K  - 1) 


,(2K  - 1)(K  + 1) 


&)}’] 


(90) 


B 


a2P 

ah* 


« !£  [4<KJ  - . (Ka  . l)u  . 2k  |^] 

Kua  L(2k  - i)uf  Vat/  i atah  at  J 


(91) 


with  B * (7K  - 5)(K  + l)8u2/4(K  - 1),  for  the  initial  values  of  the 
second  derivatives  of  particle  velocity  and  pressure  with  respect  to 
Lagrange  distance.  When  the  reaction  rate  In  polytropic  explosive  is 
known  as  a function  of  state  variables,  Eq.  (91)  can  be  used  to 


determine  whether  the  wave  produced  by  a constant  velocity  piston 
develops  Initially  with  a positive  or  a negative  pressure  gradient. 

The  validity  of  Eq.  (91)  was  first  tested  by  checking  that  dSp/dhS  = 0 
when  the  energy  release  rate  is  defined  by  Eq.  (67)  and  the  flow  develops 
as  a step  shock  in  pressure.  The  expression  for  d p/dh  as  calculated 
with 

(f-)”  <92> 


ax 

at 


P V 


(2K  - 1)  ^ 
(K  - 1)  Qtt 


- 2 3 

showed  indeed  that  a p/dh  * 0 when  n * (K  - 2)/(K  - 1),  and  also 
led  to  the  conclusion  that  aap/aha  < 0 when  n > (K  - 2)/(K  - 1)  and 
that  aap/3h2  > 0 when  n < (K  - 2)/(K  - 1). 

Expressions  for  d p/dh  were  also  derived  using  the  following 
reaction  rate  laws: 


ax 

at 

ax 

at 

ax 

at 


simulate  the  initiation  behavior  of  condensed  explosives.  The  form  of 
Eq.  (91)  for  the  Arrhenius  rate  law  will  not  be  given  here  because  the 
initiation  process  does  not  depend  on  the  shape  of  the  pressure  profile. 
That  this  is  so  follows  from  the  equations  governing  the  temperature  T. 
Since  the  volume  is  constant  along  the  shock  path,  the  rates  of  change 
of  temperature  and  pressure  at  the  shock  front  are  related  by  the  equation 


Aie~Ta/T(l  - X) 

(93) 

A^-^a  - X)“ 

(94) 

Aa(l  - X)n(p  - p )“ 

, C 

(95) 

the  well-known  Arrhenius  expression  for  a unimolecular 
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s«  (94)  and  (95)  were  used  by  Harper  and  Bernier  to 

(I) 


dT, 


H 


v,X 


dt 


dt 


(96) 


The  rates  of  change  of  temperature  and  pressure  along  a particle  path, 
however,  are  related  by  the  equation 
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lE  + p^JLil 

dt  C . 

P.A 


(97) 


where  C and  C . denote  the  specific  heats  at  constant  volume  and 
pressure^’  Initial lp* therefore  when  dp^/dt  = dp/dt,  dT/dt  > dTjj/dt  and 
dT/dh  < 0,  the  temperature  will  increase  faster  along  the  piston  path 
than  along  the  shock  path.  As  a consequence,  the  reaction  rate  will  be 
highest  at  the  piston  for  both  a positive  and  a negative  pressure  gra- 
dient. The  criticality  condition  for  initiation  is  then  that  the  piston 
motion  be  maintained  for  a time  longer  than  the  characteristic  reaction 
time 
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(98) 


when  T /T  » 1 and  T denotes  the  initial  shock  temperature, 
a i i 

The  situation  is  different  for  a pressure-dependent  energy  release 
rate  because  then  the  mechanism  of  initiation  depends  on  the  initial 
shape  of  the  pressure  pulse.  Consider  first  the  case  when  the  initial 
value  of  dap/dha  < 0.  The  pressure  and  reaction  rate  will  increase 
faster  along  the  piston  path  than  along  the  shock  path,  and  the  wave 
will  build  up  initially  from  behind.  But  when  dap/dha  > 0,  the 
pressure  and  reaction  rate  will  increase  faster  along  the  shock  path 
than  along  a particle  path,  and  the  wave  will  build  up  at  the  front. 

Since  dp/dh  >0  in  the  steady  wave,  however,  the  initiation  of  deto- 
nation will  exhibit  both  types  of  flow  when  the  initial  value  of  dap/dha  < 
Buildup  to  detonation  in  high  density  PETN  exhibiting  dp/dh  < 0 and 
dp/dh  > 0 was  recently  observed  by  Wackerle  and  Johnson. ^ As  the 
initial  pressure  is  increased,  the  mechanism  of  initiation  will  change 
from  one  involving  a negative  and  positive  pressure  gradient  to  one 
involving  only  a positive  pressure  gradient.  To  be  more  specific,  con- 
sider the  cases  when  the  energy  release  rate  is  governed  by  Eqs.  (94) 
and  (95),  and  Eq.  (91)  becomes 


and 
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(100) 


Equation  (99)  shows  that  dap/dha  < 0 when  p /p  » 1,  When  dA/dt 
is  small,  the  energy  release  rate  is  highest  at  the  piston,  and  the 
criticality  condition  for  initiation  is  that  the  piston  motion  be  main- 
tained for  a time  longer  than  the  characteristic  reaction  time 
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(101) 


Equation  (10G)  shows  that  dap/dha  > 0 when  the  initial  pressure  is 
less  than  p^  defined  by  the  condition  1 - p /p^  = (K  + 1)/Kn.  When 
this  condition  is  satisfied,  the  wave  will  buifd  up  to  detonation  with 
a positive  pressure  gradient  as  shown  by  the  results  of  Bernier's  numeri- 
cal calculations****  This  type  of  buildup  should  also  be  observed  in  the 

former  case  when  p /p.  ~ 1. 

a i 

j 2 

Calculations  of  the  initial  values  of  d p/dh  at  various  shock 
velocities  and  pressures  were  performed  using  pressure-dependent  energy 
release  rates;  a general  equation  of  state  in  which  the  unreacted 
explosive  has  a Hugoniot  curve  of  the  form  U = a + bu^  and  the  detona- 
tion product  gases  have  a poly tropic  equation  of  state;  and  the  detonation 
properties  of  PBX  9404-03,  case  Composition  B-3,  and  cast  TOT.  For  all 
three  explosives  the  initial  values  of  dBp/dha  are  negative  at  low 
initial  shock  velocities  and  then  become  positive  at  higher  shock  velocl- 
ties.  The  change  in  sign  of  d p/dh  occurs  at  a shock  velocity  sig- 
nificantly below  the  detonation  velocity  for  PBX  9404-03,  and  at  a shock 
velocity  significantly  higher  than  the  detonation  velocity  for  cast  TOT, 
The  shock  velocity  at  which  the  initial  value  of  dap/dha  becomes  posi- 
tive thus  increases  as  the  sensitivity  of  the  solid  explosive  to  shock 
initiation  decreases.  Since  d p/dh  remains  negative  for  initial 
shock  velocities  exceeding  the  steady  state  detonation  velocity  in 
cast  TOT,  the  calculations  predict  the  formation  of  a reactive  pressure 
pulse  at  the  piston  face;  this  pulse  eventually  overtakes  the  initial 
shock  front  and  causes  detonation.  This  type  of  initiation  has  been 
observed  in  homogeneous  liquid  explosives  and  in  high  density  PETN  and 
XTX-8003.  The  buildup  to  detonation  in  PBX  9404-03  and  Composition  B-3 
is  typical  of  heterogeneous  solid  explosives  in  which  the  initial  shock 
front  is  strengthened  and  accelerated  without  the  appearance  of  a 
velocity  overshoot. 

Smooth  buildup  to  detonation  with  (dap/dh2)i  > 0,  resulting  from 
the  impact  of  a flying  plate,  will  now  be  discussed.  For  convenience, 

3/ 


25 


the  flow  is  assumed  to  satisfy  the  condition 


d*p/dhs  > 0 


so  that 


3 u/dhdt  < 0,  and  the  pressure  on  the  piston  is  assumed  to  attain  a 
maximum  during  the  course  of  the  reaction.  The  relationship  between  the 
particle  velocity  gradient  and  the  energy  release  rate  when  the  pressure 
on  the  piston  is  a maximum  is  obtained  by  setting  dp/dt  =0  in  the 
energy  equation,  rewritten  as 


Sr  dp 
dh  St 


ca  du  r SX 

v dh  + Vo  **  dt 


(102) 


The  particle  velocity  gradient  on  the  piston  is  also  assumed  to  be 
positive  during  the  course  of  the  reaction.  The  curves  of  constant 
pressure  were  sketched  to  gain  a better  understanding  of  this  type  of 
flow,  and  a critical  energy  was  defined  for  the  initiation  of  detonation. 


Since  d2p/dha  >0  at  the  piston,  the  pressure  increases  faster 
along  the  shock  path  than  along  the  piston  path,  and  the  flow  develops 
with  a positive  pressure  gradient.  The  curves  of  constant  pressure  in 
the  (t-h)  plane  for  such  a flow  are  shown  schematically  in  Figure  2 
where  OCS  represents  the  shock  path,  OOP  represents  the  piston  path, 
and  DC  represents  the  locus  of  pressure  peaks. 


Particles  that  enter  the  wave  before  t attain  their  maximum 
pressure  along  DC,  but  particles  that  enter  ?he  wave  after  t attain 
their  maximum  pressure  at  the  shock  front.  Qualitative  features  of  the 
curves  of  constant  pressure  follow  directly  from  the  identity 
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where  C * (dh/dt)  denotes 
Ve  are  interested  ii  the  case 
sign  of  C is  determined  by 
a particlePpath,  and  Cp  < 0 
dp/dt  < 0.  When  dp/dh  « 0, 
and 


the  slope  of  a curve  of  constant  pressure. 

when  dp/dh  J 0.  When  dp/dh  > 0,  the 

the  time  rate  of  change  of  pressure  along 

where  dp/dt  > 0,  but  C >0  where 

however,  (dt/dh)  * 0 wEere  (dp/dt)  f 0, 
P 

(104) 


at  a singular  point  where  dp/dt  * 0.  We  are  now  in  a position  to  con- 
sider the  curves  of  constant  pressure  emanating  from  the  piston.  Since 
dp/dh  « o at  the  piston,  the  curves  of  constant  pressure  Intersect  the 


piston  path  at  right  angles  except  at  the  singular  point  D where 
Sp/St  = 0.  The  curves  emanating  from  the  piston  before  t^  have  nega- 
tive slopes  and  intersect  the  shock  path  because  the  pressure  is  higher 
at  the  shock  than  at  the  piston.  But  the  curves  emanating  from  the 
piston  after  t^  have  positive  slopes  and  do  not  intersect  the  shock 
path  because  the  pressure  is  decreasing  along  a particle  path  while  it 
is  increasing  along  the  shock  path.  When  S p/St  <0  at  D,  D is  a 
double  point,  and  the  slopes  of  the  curves  passing  through  it  are 
determined  by  Eq.  (104). 


FIGURE  2 


SCHEMATIC  DIAGRAM  OF  THE  LINES  OF  CONSTANT  PRESSURE  FOR 
THE  INITIATION  OF  DETONATION  UNITS  Op/dh)  > 0 


Consider  now  the  locus  of  peaks  DC  where  dp/dt  = 0 and  Dp/Dt  > 0 


Since 

Dp  dp  Dh 
Dt  “ 9h  Dt 


(105) 


along  DC,  the  slope  of  DC  must  be  positive  except  at  the  double  point, 
where  it  must  be  zero.  Moreover  since  dp/dh  >0,  it  follows  that  C = 0 
along  CD  for  t > t^  and  consequently  the  curves  of  constant  pressure 
intersecting  CD  are  tangent  to  the  particle  paths  as  shown  in  Figure  2. 

The  curves  of  constant  pressure  emanating  from  the  shock  path  after 

t will  now  be  considered.  It  is  evident  that  these  curves  cannot  inter- 
c 

sect  the  shock  path  again  because  the  pressure  increases  along  the  shock 
path  as  the  wave  develops  into  a detonation.  Particular  attention  should 
be  given  to  the  curve  with  the  Chapman-Jouguet  (CJ)  pressure  shown 
schematically  as  JJ ' in  Figure  2.  JJ  must  become  parallel  to  the 
shock  front  and  coincide  with  the  CJ  characteristic  and  the  line  of 
complete  reaction  as  the  reaction  zone  approaches  a steady  state.  It 
consequently  separates  the  steady  flow  from  the  unsteady  flow,  and  all 
curves  of  constant  pressure  emanating  from  the  unsteady  shock  above  J 
must  become  parallel  to  the  shock  front. 

A critical  energy  E for  the  shock  initiation  of  an  explosive 
can  be  defined  as 

E = i pudt  (106) 

c J o 


where  p and  u are  the  pressure  and  particle  velocity  along  the  piston, 
the  half  is  included  to  account  for  the  conversion  of  work  into  kinetic 
and  internal  energy,  and  the  critical  time,  Of,  denotes  the  time  of 
formation  of  the  forward- facing  characteristic  that  eventually  becomes 


the  CJ  characteristic  in  the  steady  detonation  wave.  Note  that  the 


critical  energy  is  defined  as  a change  in  the  internal  energy  rather 
than  as  the  change  in  total  energy  specified  by  Walker  and  Wasley.15 
Since  former  work  on  initiation1  suggests  that  this  critical  charac- 
teristic lies  in  a centered  reactive  rarefaction  fan,  Of  can  in 
principle  be  determined  by  finding  the  position  on  the  piston  path  where 
the  second  derivative  of  the  characteristic  becomes  zero.  Let  s 
denote  the  slope  of  a forward-facing  characteristic,  so  that 


(dr/dh) 


The  condition  determining  0!  can  then  be  written  as 


+ /ScA 

dt  “ Vo(  at 


0 = °* 


(108) 


Attempts  were  made  to  calculate  p,  u,  and  On  along  the  piston 
to  determine  Ec  for  a poly tropic  medium  with  a simple  reaction  rate 
law  but,  since  approximations  must  be  made  unless  a complete  solution 
to  the  reactive  flow  problem  i»  knovrn,  a simpler  approach  to  the  critical 
energy  problem  was  undertaken. 


* The  alternative  approach  was  to  apply  Eq.  (106)  to  the  steady  state 
wave  and  assume  that  the  amount  of  energy  required  to  initiate  a steady 
detonation  with  no  delay  time17  is  sufficient  to  initiate  the  explosive 
subjected  to  other  initial  conditions.  Another  way  to  look  at  this 
assumption  is  to  say  that,  if  there  is  a critical  energy  for  initiation, 
its  value  will  be  determined  by  the  properties  of  the  steady  state  wave. 
Since  the  centered  rarefaction  fan  in  a steady  wave  emanates  from  the 
CJ  point,  the  critical  time  a in  the  steady  wave  is  taken  as  the  reac- 
tion time.  Evaluating  Eq.  (106)  for  a steady  state  detonation  wave  leads 
to  the  following  equation  for  the  critical  energy 
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2 <P0/P> 


(109) 


where  D denotes  the  detonation  velocity,  Z denotes  the  reaction  zone 
length,  and  the  density  terms  in  brackets  denote  values  averaged  over 
the  reaction  zone.  Equation  (109)  for  the  critical  energy  takes  into 
account  the  initial  state  and  the  equation  of  state  of  the  unreacted 
explosive,  the  equation  of  state  of  the  reactive  mixture,  and  the  kinetics 
and  thermochemistry  of  the  exothermic  reaction  supporting  detonation. 
However,  it  gives  values  larger  than  those  determined  experimentally  by 
Walker  and  Wasley.  This  disparity  suggests  that  calculations  of  critical 
energy  with  the  steady  state  wave  must  account  for  the  fact  that  the 
completion  of  reaction  on  the  rear  boundary  is  not  a necessary  condition 
for  initiation  of  detonation. 
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VI  THE  REACTIVE  SHOCK  PROBLEM 


Work  was  continued  on  the  reactive  shock  problem  with  a prescribed 
energy  release  rate.  Solutions  for  buildup  to  detonation  were  sought  to 
provide  a means  of  demonstrating  the  dependence  of  initiation  on  the 
energy  release  rate  and  on  the  rear-boundary  conditions,  and  a means  of 
calculating  the  critical  energy  for  initiation  defined  in  Section  V. 

The  new  methods  tried  for.  constructing  solutions  for  reactive  shocks 
were  unsuccessful  however. 


It  is  convenient  to  write  the  flow  equations  in  terms  of  the 
Lagrange  time  T introduced  earlier  in  this  report  as 


dv  _ vo  9u 
5t  " U St 

du  _ vp  dp 

St  " 'u  St 

dp  _ /c~\  dv  r ax  . 

dt  \v/  dt  + v ^ dt 


(110) 


(111) 


(112) 


Our  problem  for  a prescribed  energy  release  rate  is  to  find  the  volume 
field  v(t,T),  the  particle  velocity  field  u = u(t,T),  and  the  pres- 
sure field  p(t,T),  that  satisfy  the  differential  equations  (110)-(112) 
along  particle  paths  and  satisfy  the  Rankine-Hugoniot  jump  conditions 
along  the  shock  path  where  t * T,  Previous  attempts  to  solve  this 
problem,4  with  the  explosive  described  by  Eq.  (23),  and  the  energy 
equation  written  as 


dr  _dp 
dT  dt 


W1* 


(113) 


made  use  of  the  fact  that  the  v(t,T),  u(t,T),  and  p(t,T)  fields  can 

be  obtained  from  either  the  particle  velocity  gradient  or  the  pressure 
gradient  by  integrating  Eq.  (38),  Eq.  (110),  and  Eq.  (Ill)  when  the  energy 
release  rate  at  the  shock  front  is  prescribed.  Consider,  for  example, 
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the  case  when  the  form  of  du/dT  Is  assumed.  Integration  of  Eq.  (38) 
gives  the  shock  pressure,  the  shock  particle  velocity  follows  from 
Eq.  (37),  and  the  shock  path  can  then  be  obtained  by  Integrating  Eq.  (35). 
Integration  of  Eq.  (10)  along  a particle  path  from  the  shock  front  gives 
the  volume  field,  and  Integration  of  du/dT  and  Eq.  (Ill)  along  an 
Isochrone  from  the  shock  path  gives  the  corresponding  particle  velocity 
and  pressure  fields.  Since  the  energy  equation  was  not  used,  however, 
a solution  for  a reactive  shock  Is  obtained  only  If  the  assumed  form  of 
du/dT  is  compatible  with  the  energy  release  rate  and  the  calculated 
flow  satisfies  Eq.  (113).  A similar  argument  applies  also  to  the  case 
when  the  form  of  dp/dT  is  assumed.  The  problem  can  therefore  be 

regarded  as  that  of  finding  the  particle  velocity  gradient  that  is  com- 
patible with  the  prescribed  energy  release  rate.  This  approach  is  prac- 

tically intractable,  however,  because  the  flow  is  related  to  the  reaction 
through  the  energy  equation,  and  there  is  no  apparent  way  that  this 
equation  can  be  used  to  determine  the  particle  velocity  gradient. 


Since  buildup  to  detonation  is  of  particular  interest,  most  atten- 
tion was  given  to  the  final  attainment  of  steady  state  flow.  In  one 
approach,  the  equations  of  motion  were  transformed  so  that  their  integrals 
consist  of  steady  state  terms  and  unsteady  state  terms  that  vanish  as 
the  steady  state  is  attained.  The  independent  variables  were  changed 


from  ( t, t)  to  (t,£  = T - t)  so  that  the  partial  derivatives 
transform  as  follows 


(114) 


(115) 


and  Eqs.  (110)  and  (111)  become 


(116) 


(117) 


Formal  Integration  of  Eqs.  (116)  and  (117)  with  respect  to  £ gives  the 
equations  expressing  the  conservation  of  mass  and  momentum  for  unsteady 
flow  as 


(118) 


(119) 


The  corresponding  equation  for  the  energy  is  most  conveniently  derived 
by  integrating  the  first  law  of  thermodynamics  written  as 


de  dpv  dp  . 

H + H m v H 


(120) 


The  combination  of  Eqs.  (120),  (117),  and  (118)  leads  to  the  equation 


de  dpv  , . _du  du  dp  du  r /d u\ 

h * 3*-  * <u  - - u ar-  - v ST  - ae  Jo 
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which  can  readily  be  integrated  to  give  the  equation  expressing  the 
conservation  of  mass,  momentum,  and  energy  as 

. - eo  * iw  - p0v„  * i(u  - B)a  - Jo"  = 

Equations  (118),  (119),  and  (122)  can  be  regarded  as  perturbed  Rankine- 
Hugoniot  equations  because  thoy  reduced  to  these  equations  as  the  flow 
attains  a steady  state  and  the  derivatives  (du/dTK  = 0 and  (dp/dT)^ 
becomes  zero.  Attempts  to  u&e  these  equations  in  aSperturbation  analy- 
sis to  construct  solutions  '.'or  the  final  approach  to  the  steady  state 
were  unsuccessful. 

In  the  second  approach  to  the  reactive  shock  problem,  the  energy 
equation  was  integrated  rather  than  used  as  a compatibility  condition 
as  in  the  previous  treatment.4  Although  the  pressure  and  volume  fields 
generated  from  the  energy  equation  are  compatible  with  a prescribed 
energy  release  rate,  there  is  no  guarantee  that  they  are  associated 


with  a unique  particle  velocity  field  because  the  continuity  and  momentum 
equations  were  not  used  in  their  construction.  The  continuity  and 
momentum  equations  were  therefore  combined  to  obtain  a condition  that 
the  pressure  and  volume  must  satisfy  in  order  to  be  a solution  to  a 
reactive  shock  problem.  This  compatibility  condition  was  derived  as 


H 3p  a day 
u ——  t*-  + p rr—s 
M dT  3t  H St* 
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by  eliminating  dSu/dtdT  from  the  equations  obtained  by  differentiating 
Eq.  (110)  partially  with  respect  to  t and  differentiating  Eq.  (Ill) 
partially  with  respect  to  T. 


A simple  case  was  considered  first.  An  expression  relating  the 
volume  and  reaction  was  chosen  as 
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so  that  the  energy  equation  could  be  integrated  without  difficulty,  and 
the  reaction  was  assumed  to  have  no  activation  energy.  The  reaction 
rate  was  written  formally  as 


ax 

at 


i Ml  - X) 
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with  R(l)  * 1 and  R(0)  = 0,  so  that  aX/3t  = - aX/dT  and  the 
reaction  rate  at  the  shock  front  is  constant.  Integration  of  the 
energy  equation  gives  the  equation  for  the  pressure  field  as 


[F(A)  - F(0)] 


(126) 


4(K  - l)G(X) 


where  p^  denotes  the  spike  pressure  in  a CJ  wave  and  F is  defined 
by  the  equation  dF/dX  * G*-1.  Integration  of  the  differential  equation 
along  the  shock  path  gives  the  equation  for  the  shock  pressure  as 


g(0)t/aj  + P 


(127) 


HO 


Equations  (126)  and  (127)  then  give  the  equations  for  the  pressure  field 
as 


PD i.  VT>  ' PD 

p s ~ (l  + (l  - X)  ) + ~ 
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The  fact  that  Eqa.  (124),  (131),  and  (132)  do  not  satisfy  Eq.  (123) 
leads  to  the  conclusion  that  Eq.  (124)  cannot  be  satisfied  In  a tine- 
dependent  wave  that  builds  up  to  a steady  state.  In  other  words, 

Eq.  (124)  Is  too  restrictive  and  nust  be  generalized  to  obtain  a more 
realistic  treatment  of  unsteady  flow. 


Other  attempts  to  construct  solutions  for  the  final  stages  of 
initiation  were  based  on  the  assumption  that  the  fall  in  pressure  along 
a particle  path  is  proportional  to  the  fall  in  pressure  along  a particle 
path  in  the  steady  state  wave.  The  equation  expressing  this  condition 


is  obtained  from  the  equation  for  the  pressure  in  a CJ  wave 
P * pd/2(1  + (l  - A*)  as 


(133) 


with  E an  arbitrary  function  of  T.  A unimolecular  reaction  with  no 
activation  energy  was  considered  to  simplify  the  integration  of  the 
energy  equation.  In  this  case,  Eq.  (133)  reduces  to  the  equation 
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since  1 - X = e^T  Integration  of  Eq.  (134)  gives  the  equation 

for  the  pressure  field,  which  is  written  for  convenience  as 


B + E 


with  B another  arbitrary  function  of  T.  Integration  of  the  energy 
equation  gives  the  corresponding  volume  field  as 

*' = TFTuI^  - E *“  ) ‘ U36) 


(135) 


Evaluating  Eq.  (136)  at  the  shock  front,  however,  where  v = v = 
[(K  - 1)/(K  + l)]v0  and  t ■ T gives  the  following  relationship 
between  E and  B, 


KB  = (K  - 1)E2  + E . 


(137) 


The  equations  for  the  pressure  and  volume  can  therefore  be  written 
in  terms  of  the  arbitrary  function  E as 


P = 


(K  - 1)E‘ 


T-ty 

Jl  20! 

+ lit  + e 


(138) 


and 


v = 


vo. 


(K  + 1) 


((C  - 1)  + - 1 - e 

E 


2a 


(139) 


The  equation  for  the  shock  pressure  follows  from  (138)  as 


rH 


P 

-£  ((K  - 1)E2  + (K  + 1)e) 


(140) 


and  the  differential  equation  for  shock  pressure  leads  to  the  following 


equation  for 


E 


* 

; 


dE 

dt 


Cl  - E) 

6a(l  + [(K  - 1)/(K  + 1)]e) 


(141) 


The  question  now  arises  whether  Eqs.  (138)  and  (139)  satisfy  the  com- 
patibility condition.  The  equations  obtained  by  differentiating  Eqs. 
(138)  and  (139),  together  with  Eq.  (141),  show  that  Eq.  (123)  is  not 
satisfied.  Here  again  the  assumption  made  about  the  flow  was  not 
realistic  enough  and  we  were  unable  to  find  a solution  for  the  final 
buildup  to  detonation. 


VII  RESULTS  AND  RECOMMENDATIONS 


The  results  of  the  theoretical  study  of  reactive  shock  waves 
presented  in  this  report  lead  to  a better  understanding  of  the  initiation 
of  detonation  in  condensed  explosives. 

Various  aspects  of  shock  initiation  were  considered.  The  differ- 
ential equation  governing  a shock  discontinuity  was  used  to  determine 
different  conditions  for  a single  shock  trajectory  for  buildup  to 
detonation.  One  of  these  conditions  for  single  shock  buildup  produced 
by  a flying  plate  was  used  to  construct  a solution  for  the  type  of  flow 
observed  by  Kennedy10  during  the  early  stages  of  initiation  in  PBZ  9404. 

Other  properties  of  initiation  induced  by  a constant  velocity  piston 
or  a flying  plate  were  considered.  The  equations  relating  the  initial 
flow  to  the  initial  energy  release  rate  in  such  a wave  were  derived. 

These  equations  provide  a means  of  determining  the  variation  of  energy 
release  rate  with  pressure  at  the  front  of  the  wave  when  the  Hugoniot 
curve  of  the  explosive  is  known.  Conditions  were  also  determined  for  the 
shock  produced  by  a constant  velocity  piston  to  accelerate  with  either 
a positive  or  a negative  pressure  gradient.  These  conditions  are  important 
because  they  demonstrate  how  the  mechanism  of  initiation  depends  on  the 
energy  release  rate,  the  sound  speed,  and  the  relationship  between  these 
quantities.  A critical  energy  for  initiation  was  defined  for  the  case 
when  the  wave  builds  up  with  a positive  pressure  gradient.  Vork  on  the 
reactive  shock  problem  was  continued,  and  integral  relationships  for 
unsteady  flow  were  derived  as  generalized  Rankine-Hugoniot  equations 
without  making  approximations.  Attempts  to  construct  explicit  solu- 
tions for  the  initiation  of  steady  state  waves  so  that  critical  energies 
could  be  calculated  were  however  unsuccessful. 

Steady  state  detonation  was  considered,  and  a perturbation  analysis 
for  a simple  reaction  was  carried  out  to  determine  a sufficient  condition 
for  its  stability. 

Additional  studies  that  would  extend  the  work  presented  in  this 
report  are  recommended  to  Improve  our  present  understanding  of  the 
initiation  of  detonation  in  condensed  explosives.  A general  criticality 
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Appendix 


STABILITY  OF  THE  STEADY  DETONATION 


i 

l 


An  objective  of  the  research  reported  here  has  been  to  separate 
initial  conditions  leading  to  steady  detonation  from  those  producing 
extinction  or  some  sort  of  irregular  burning.  For  initial  conditions 
close  to  the  steady  state , the  problem  is  one  of  finding  the  stability 
boundary  where  perturbations  in  parameter  values  on  one  side  of  the 
boundary  lead  to  steady  detonation  while  those  on  the  other  side  do  not. 


A first  attempt  at  locating  the  stability  boundary  is  made  below 
by  examining  whether  a small  initial  irregularity  will  grow.  The 
irregularity  at  a point  in  the  unbumt  material  is  swept  over  by  the 
shock.  The  perturbation  wave  in  the  reacting  flow  behind  the  shock  is 
headed  by  a sound  wave  or  characteristic.  It  is  found  that  the  effect 
of  the  initial  disturbance  dies  away  along  the  characteristic  if  K , 
the  polytropic  index,  is  small  enough;  the  critical  K depends 
on  n , the  order  of  the  reaction.  If  K is  larger  than  this  value, 
the  initial  disturbance  does  increase  with  time,  at  least  part  way  along 
the  characteristic,  so  that  an  instability  can  arise.  The  results  is 
sketched  for  a particular  case  in  Figure  A-l. 


Steady  State  ^ 

The  gas  dynamic  equations  for  one-dimensional,  adiabatic,  reactive 
flow  read 


> + p"» 
Pu  + p_ 


0 

0 


p + KpUg  » (K  - 1)QPR 


.-.r 

0m 


X . 


R 


- - - ■ 

- xji  : . . 

- ’4-"  ^ 

, 

- . 

v 

< r ‘ 

•:vP-r^aBgf- 

. h 


8J®| 


where  p is  the  gas  density,  p the  pressure,  u the  particle  velocity, 
X the  fraction  of  gas  reacted  at  rate  R,  K the  polytropic 
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index*  and  Q the  heat  released.  The  Eulerian  coordinates  are  distance 
x and  time  t with  the  superior  dot  denoting  the  material  time 
derivative. 


* 


\ 


I 


Solutions  of  the  equations  representing  steady  states  sure  those 
where  P , p , u , and  X are  functions  of  (t  - x/D)  only,  where 
D is  a constant  speed.  If  the  steady  state  is  headed  by  a strong  shock 
moving  at  speed  D into  quiescent  gas  of  density  Po  , the  jump  conditions 
give  just  behind  the  shock 

u * 2D/(K  + 1) 

p « Po  <K  + 1)/(K  - 1) 

p * PDU 

X * o 

♦ 

The  corresponding  steady  state  is 

P « PoD/(D  - u) 
p = PoDU 

X a u£D  - (1C  + 1)  u/2]/(IC  - 1)Q 

The  Chapman- Jougue t condition,  that  D * u + c where  the  reaction  is 
complete  (X  = 1)  , is  satisfied  if 

Da  » 2((CS  - 1)Q 


Here  c is  the  sound  speed  given  by 

Ca  a Ku(D  - u) 

The  above  steady  state  equations  relate  the  dependent  variables  to 
one  another.  Their  relation  to  the  coordinate  locations  is  found  from 
the  reaction  equation 

sas-  V * 

This  equation  can  be  expanded  to  read 

2(*  + 1)(D  - u)  [D  - (K  + 1)u3uVd*  “ H 


S-t: 
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where  the  prime  denotes  differentiation  with  respect  to  (t 
Integration  gives 

u 

t - x/D  - J 2(IC  + 1)(D  - u)  [D  - (K  + l)u]du/D*R 


- x/D). 


The  particular  form  for  R used  below  is 
r =.  p"“l(i  - X)?a 

where  OL  and  n are  constants.  This  form  approximates  the  rate 
expression  for  a reaction  of  order  n with  simple  pressure  dependence. 
Since 

1 - X * [D  - (K  + l)u]a/D8 

one  sees  that  this  reaction  is  complete  in  a finite  time  only  if  n < 1. 


Perturbation  Equations 

For  a small  perturbation  of  the  steady  state, 

P * P + ep,  u « u + Cu,  etc., 

where  the  bar  indicates  the  steady  state  function  given  above,  the  tilde 
denotes  a quantity  of  ordinary  magnitude,  and  e is  a small  constant. 
Substitution  in  the  gas  dynamic  equations  gives,  to  the  first  order  in 


+ up  + up  + K(pu  + pu  ) ■ (K  - 1)  Q(PR  + PR) 

t X X X X 


where  the  tildes  have  been  omitted  from  the  perturbation  quantities  for 
simplicity. 


Along  particle  paths,  linear  combinations  of  the  equations  in 
characteristic  form  read: 

A + uA  + uA  = R 

t X X 

- ca  (p^  + up  ) + p^  + up  - ca  (up  + Pu  ) + p U + tCpu 
t x t X X X X X 

a (K  - 1)Q  (pR  + pR) 

along  particle  paths,  and 

p + (u  ± c)p  ± cP[u  + (u  ± c)u  ] 

X XX  X 

± c(p[u  + uu  ] + Pu  u)  + p u + Kpu 

t X X X X 

- <K  - DQ  (PR  + pR) 


along  forward  (+)  and  backward  (-)  sound  waves. 


The  perturbation  at  the  shock  is  found  in  the  same  way: 
u » 2D/(K  + i) 

P a 0 
p"  a 2pDu 
A = 0 


Irregularity  at  a Point 


A small  irregularity  is  supposed  present  at  a point  in  the  unburnt 
material.  When  the  shock  reaches  this  point  it  will  experience  a 
velocity  perturbation  SD.  Corresponding  perturbations  u,  p,  p,  and 
A arise  according  to  the  shock  equations  given  above.  The  expressions 
are  made  precise  by  setting  Dal  and  putting  the  magnitude  of  the 
perturbation  into  6.  The  only  requirement  is  that  e is  small  enough 
that  it  is  reasonable  to  neglect  e with  respect  to  1. 


The  first  effect  of  the  irregularity  is  a discontinuity  propagated 


along  the  backward  running  sound  wave, 
sound  wave,  the  equation 


Sihce  particle  paths  cross  this 


X + uX  + uX  « R 

t X X 

holds  across  the  wave  in  spite  of  the  jump  in  the  value  of  u and  per- 
haps of  R.  One  sees  that  the  equation  remains  satisfied  by  an  abrupt 
change  in  the  cross-derivative  of  X, 


if  needed,  and  not  by  any  change  in  X itself.  The  perturbation  X 
therefore  retains  just  behind  the  sound  wave  the  value  zero  that  it  had 
in  the  steady  flow  ahead  of  the  wave.  On  the  wave 


X = 0 


The  other  equation  valid  along  particle  paths  may  be  written 
(p  - cs  p)  + u(p  - ca  P)  + P(ca  + uca  ) - ca  (up  + pu  ) + p u + Kpu 

t X t X X x x 

= (K  - 1)Q(PR  + PR) 

By  the  same  argument  as  above,  one  finds  that 

P - C2  P SB  0 

on  the  sound  wave. 


Frcxn  the  forward  running  sound  wave  in  the  same  way, 


The  above  relations  valid  along  the  backward  running  sound  wave 
permit  the  elimination  of  p,  p,  and  X from  the  equation  for  the 
wave.  One  Obtains 


u + (u  - c)u  + (u  /4KcDu)Mu  * 0 


where 


M SB  KauD  + C[(K  - 2)5  + 2(1  - Ka)uJ  + CK53«/(K  + l)u  3(R/u) 


Since  u has  the  same  sign  as  R , which  is  positive,  the 

..  , . ..  .. 

coefficient  of  Mu  in  the  differential  equation  is  positive,  (hie  sees 
that  the  perturbation  at  the  shock  grows  with  time  along  the  sound  wave 
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if  M is  negative  and  decays  if  M is  positive.  Since  M varies 
with  u , decay  and  therefore  stability  are  assured  only  if  M is 

positive  for  the  whole  range  of  values  of  u.  This  stability  criterion 
is  examined  for  a particular  rate  expression  below. 

Stability  of  a Particular  Reaction 


For  the  reaction 


R = Pn“1(l  - X)n/« 


one  finds 


R/R  = (n  - 1)  (p/p)  - nX/(l  - X) 


Since  along  the  backward  running  sound  wave  propagated  from  the  initial 
disturbance 


X = 0 


one  finds 


(P/c)u 


R/u  = (n  - l)R/c 


Use  of  the  expressions  for  these  steady  state  quantities  in  terms  of  u 
gives  finally 

M = KauD  + c[(K  + n - 3)5  + (K  + l)u(3  - n - 2*)] 

An  investigation  of  the  equation  of  M for  u on  its  range  D/(K  + 1)  ^ 
u £ 2D/(k  + 1)  is  now  required. 

If  one  sets  M » 0 , one  finds  that  there  are_generally  three  roots 
for  u.  One  of  these  roots  is  at  the  C-J  value,  u * D/(K  + 1),  for 
any  values  of  K and  n.  There  is  a second  root  at  u = D/(K  + 1)  if 

n = 3 + K(K  - 3) 


There  is  a root  at  the  shock,  u = 2D/(K  + 1)  , if 


n - 3 + k{-2  ± [<K  + 1)/(K  - 1)]*} 
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